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Abstract 

: We present analytic and numerical results for the evolution of currents on super- 
conducting strings in the classical U(l) X U (1) model. We derive an energy functional 
for the currents and charges on these strings, establishing rigorously that minima 
should exist in this model for loops of finite size (vortons) if both charge and current 
are present on the worldsheet. We then study the stability of the currents on these 
strings, and we find an analytic criterion for the onset of instability (in the neutral 
limit). This limit specifies a lower maximal current than previous heuristic estimates. 
We conclude with a discussion of the evolution of loops towards their final vorton 
state in the model under consideration. 
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1 Introduction 



Topological defects are a class of exact solutions in field theories whose stabil- 
ity is enforced by topological reasons. In particular, strings, the class of defects 
associated with a non-trivial first homotopy group of the vacuum manifold, 
have been widely studied, since they seem to appear in a variety of gener- 
alisations of the Standard Model (GUTs, SUSY etc.). They are for example 
associated with the spontaneous symmetry breakdown of a U(l) symmetry, 
like the axion or baryon symmetry (global), or electromagnetism (local) in 
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superconductivity. Strings could therefore appear during a cosmological phase 
transition, and are prime candidates for a number of astrophysical puzzles, 
like the dark matter of the universe or the origin of the most energetic cosmic 
rays (for a review of cosmic defects, see ref. [1]). 

In this paper, we shall be interested in a class of defects where the string's 
field is coupled to another scalar field, which allows the build up of charge and 
currents on the worldsheet. Because of the non-dissipative properties of the 
currents on these strings, they are called superconducting. 

The plan of this paper is as follows: in the next section, we will discuss the 
Lagrangian under consideration, and show how the amount of charge and 
currents is limited; in the following part, we will carry out an analysis of the 
energy functional of the condensate, to get a very simple analytic expression 
for it. We will then use our results to discuss the possibility of forming stable 
loops of superconducting strings. Finally, we will study the stability of currents 
on the string's worldsheet and derive an exact analytic result for the onset of 
instability, which we shall illustrate by numerical simulations obtained from 
our full 3D field theory code. 



2 The field theory model 



The model under consideration here is the original U(l) x U(l) model, first 
proposed by Witten [2], and based on the Lagrangian 



-^(w\ 2 - vir - m 2 w\ 2 ■ a) 

where (f> and a are complex scalar fields and A^, A CT , r]^, rj a and j3 are posi- 
tive constants. We can arrange the parameters in this model such that the 
ground state has the ^-symmetry broken (|0| = r]^ ^ 0), while a remains 
symmetric (\a\ = 0). Under these circumstances, vortex solutions exist in the 
0-field. Here, we shall assume we have a vortex-string lying along the z-axis; 
in azimuthal coordinates, this solution takes the form, 

</>(r,6,z) = \<l>\(r)e», (2) 

with |0|(O) = at the vortex centre and |0| — > rj^ as r — > oo. 

We wish to consider the conditions under which a condensate in the cx-field 
can emerge in the core of the string. This condensate can also carry currents 
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and charges along the string, so we will represent it by the following ansatz 
which describes the dependence of these excitations on z and t: 



a= \a\(r) exp^ t+kz \ 



(3) 



We can see that the presence of charge and current causes a change in the 



Lagrangian SC 



k )\a\ , which alters the mass term [3,4], and the vac- 



uum expectation value of a; indeed, we can rewrite the Mexican hat potential 
for Icrl as: 



id 2 -k+£^ 2 -* 2 f 2 



(4) 



From the expression (4), we see that the 0-symmetry will remain broken in 
the vacuum, with \(f>\ = rj^, as long as we satisfy 



(5) 



If we want the cr-symmetry to be remain unbroken in the ground state (|cr| = 
0) , then we need to ensure that the mass term for a is positive at infinity, that 
is, 

< = (3vl-lWa-^ 2 + k 2 >0. (6) 

Now at the centre of the string we have \<f>\ = 0, so this mass term can become 
negative and a can develop a non-vanishing expectation value, provided that 



-Kr,l + cu 2 - k 2 > . 



(7) 



However, the condition (7) is not sufficient to obtain a cr-condensate within the 
string, since we also have to consider its gradient energy cost. To determine 
this we follow the analysis of Haws et al. [5] and study the stability of the 
trivial solution a — with perturbations of the form a = \a(x,y)\e lut , with 
Id << 1. The field equation, using the modified Mexican Hat potential (4), 
then becomes 

-drSTa - + 2(UJ \~ k2) ] ~ P^r 2 ) a = u 2 a , (8) 

which is an eigenvalue equation with a harmonic oscillator-like potential. 
Hence, we know the value of the ground state eigenvalue vq, and we see that 
the trivial solution will be unstable for Vq < or, equivalently, 

k 2 -u 2 K^Kri-JpX^l. (9) 
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Figure 1. Profiles of the vortex and the condensate fields, with no charge and current 
(solid lines). The effects of adding a significant current (dashed) or charge density 
(dotted) is also illustrated; the (anti-) quenching effect of current (charge) is clear. 



If (5), (6), and (9) are satisfied, a will form a condensate of width S a ~ m~ x , 
with two conserved quantities, the usual Noether charge Q and a topological 
charge N, where: 

Q = u J dz J dS\a\ 2 , N = J dz ^ , (10) 

and a current flowing along the z direction: 



J z — k 



JdzJdS\a\ 2 . (11) 



The field configuration for a vortex with a condensate can be seen in fig. 1, 
with the influence of k and uj also shown. 

To conclude this section, we note that in the bare case (uj — k — 0), (5), (6), 
and (9) lead to the following inequalities on the parameters of the theory: 

P<^, H<vl- (12) 

These equations will prove to be useful in the rest of the paper. 



3 The energy functional 



Starting from the Lagrangian (1), we can easily switch to a Hamiltonian for- 
malism. In particular, we will be interested in the following in the energy of 
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the condensate, which can be written in the form: 

s a = Jdzjds (W1 2 + (k 2 + c 2 )m 2 + {m 2 - + ^ n 4 ) • 

(13) 

This is a rather complicated expression, so in order to simplify it we consider 
the equation of motion for a: 



-d r d r a + (k 2 - u 2 )a + (/?|0| 2 - \\ c r)l)o + y M V = ■ (14) 

Multiplying (14) by the complex conjugate <r + , integrating by parts, and then 
inserting the result in (13), we find: 

£ a = -^-J d 3 r H 4 + 2 J d 3 rcu 2 \a\ 2 . (15) 

(It is easy to generalise this equation to the gauged case, making the obvious 
replacement u — > u — A t .) 

Now, if we recall the expression for the conserved charge Q (10), we can rewrite 
the energy as: 

^-^/,V H * + -g_. (16) 

This equation (16) is interesting because it expresses the energy in the conden- 
sate as a function of conserved quantities, and of the two integrated moments 
S 2 and £4, where 

£ n = / dS\a\ n . (17) 

Note that dT> n /dz = with the ansatz (3), and so we can carry out the 
^-integration in (16) for a segment of finite length L to obtain 

which is a purely analytic expression. 

This is the form of the energy functional we sought. All that is now required 
is to find an ansatz for the S n , which can be done by noticing that: 

Z n = JdS\a\ n ~ C x ^ ~ ^f, (19) 

where 5 a is the average width of the condensate, and a max is its maximum 
height in the vortex core. Now, we have the mass given by 

ml = Pvl ~ \Wa ~uj 2 + k 2 = oo 2 c -uj 2 + k 2 , (20) 
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Figure 2. Plot of the numerical values for <7 max , as a function of uj 2 



n . 



and we can see from fig. 2 that, up to a very good approximation, we can take 



a (k 2 c - k 2 + uj 2 ) 



(21) 



where k 2 ~ ^\ a T) 2 — yJpX^rj 2 , is the critical winding number density at which 
the condensate must vanish (obtained from (9)). (We note that in the gauged 
case, the picture is qualitatively the same, but that the gauge fields maintain 
the condensate against the quenching effect of the current; the height is then 
nearly constant in a somewhat wider range of uj 2 — k 2 . Ultimately, however, 
k 2 is the same for both gauged and global cases.) If we further assume that 



k c ~uj c ~ 



m , we can factorize it out in (19) to yield: 



(22) 



with v = {uj 2 — k 2 )/m 2 and S no = S n (w 2 — k 2 = 0) is the value of £„ in the 
so-called chiral case (uj = ±k). In particular, for n = 2 this reduces to: 



S 2 = E 



m 2 Q + (uj 2 - k 2 ) 



O o 

mi 



[UJ^ 



k 2 ) 



(23) 



As can be seen from fig. 3, the agreement between the numerical calculation 
and our analytic estimate (22) is remarkably good. 



4 The behaviour of superconducting loops 



Cosmic strings tend to intercommute, and therefore we can expect the copi- 
ous production of loops from a network of strings. Superconducting strings 



6 



co 2 -k 2 co 2 -k 2 

Figure 3. Plots of the numerical values for £2 and S4 (plain), and their analytic 
estimates (dashed) from (22). 

follow the same behaviour [6], but this time the charge and current due to 
the additional structure may prevent the shrinkage of the loop and stabilise 
the configuration [4]. These stable loops or vortons, effectively held up by 
their angular momentum, have been the subject of numerous studies (see, for 
example, refs. [7,8,9,10,11]). 

Using the results of the previous sections, we would like to study the behaviour 
of these loops, starting from our analytic result for the energy of the condensate 
(18). However, we must also include the energy in the vortex, which we will 
model by 8$ = fJ,L, where fj, ~ 0{rjfy is the mass per unit length of the string. 
We can then write our energy functional for the loop as 

fi^-^ + g. (24) 

We will now discuss the various possible regimes for the charges and currents 
on these loops: 

The chiral case, u 2 = k 2 : 

We will first consider this limit, where oj 2 — k 2 = (see refs. [12,13] for a study 
of chiral loops using the 2D approach of [10]). In this regime, we always have 
£ = S G and £4 = Ti 4o . Since these values remain fixed as the loop shrinks, it 
is easy to vary (24) to see that it will have a minimum for 

r2 2Q 2 2iV 2 £ 



£ o (/i-^£ 40 ) (/i-^£ 40 )' ^ 
Hence, the loop will contract (or expand) until u 2 = k 2 = (/i — iA CT £ 4o )/(2£ ). 
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For this to be well defined, we should verify that £ fe=w=0 = /i f-S 4o is 

positive. To do so, we note that this is the total energy in the bare supercon- 
ducting string, and that the profiles \4>\(r) and \cr\(r) minimise the value of 
this functional. By considering <f>(\r) and cr(Ar), differentiating the field en- 
ergy functional (13) with respect to A and setting this quantity to at A = 1, 
we obtain the result: 

£*^,=o = / d 3 r(|V0| 2 + |Va| 2 ) > . (26) 

This ensures that the loop length L 2 is indeed well defined and that chiral 
loops must have a stable minimum of their energy. (Of course, we assume 
here that chiral currents are long-lived and that the final stable loop size is 
significantly larger than the vortex width.) 



The electric and magnetic cases, uj 2 — k 2 ^ 0: 

To study these, we have to manipulate the form of the loop energy. Using the 
fact that N = kL and Q = ujY,L are conserved during the evolution of the 
loop, we can express L in the following terms: 



N 2 



UJ 



2 -k 2 \ S 2 



- 1 



(27) 



where Eqjv — Q/N. (Note that the ratio Q/N is defined at loop creation and 
is assumed to be conserved thereafter.) By substituting this back in (24), we 
get an expression for the energy as a function of uj 2 — k 2 : 



S = N 




+ 2S 



QN 



k 2 



Y 2 



£ 2 



(28) 



Here, note that E n is clearly a function of uj 2 — k 2 , since it is determined by 
the solution of the field equations, which only involve uj 2 — k 2 . 

It is also important to notice that Eqat can be expressed as: Eqjv = f E. So, 
using the terminology of Carter and Peter [10], in the electric case (uj 2 > k 2 ), 
we have S G < E < T,q N , while in the magnetic case (k 2 > uj 2 ), we have 
Tjqn < S < S . This implies that for each loop, E can only take a limited range 
of values, the bounds of which are fixed, on the one hand, by the parameters 
in the Lagrangian determining the bare E G and, on the other, by the amount 
of charge and current on the loop fixing £qjv. 



A simple analysis of (28) now tells us that, as S 2 — > T?q N (which is achieved 
for a non-vanishing value of uj 2 — k 2 ), the energy £ goes to +oo, while in the 
opposite limit as uj 2 — k 2 — > then S diverges once again, that is, provided 
we have fi — 4f X 4 > 0. Hence, as long as this condition is satisfied, then the 
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Figure 4. This figure shows the energy E as a function of lo 2 — k 2 , in the mangetic 
regime (Xqat/£ = 0.5 on the left), and in the electric regime (£qjv/X = 2 on the 
right). 

loop energy 8 must have a minimum, an equilibrium state corresponding to 
the vorton. 

As can be seen from fig. 4, the vorton state is usually reached when E ~ Sqjv, 
which means that uj/k ~ 1. This appears to justify the assertion that % — > 1 
is an attractor [4], but it does not imply that the chiral state with oj 2 — k 2 = 
is an attractor. On the contrary, a typical loop created with small charges and 
currents (\u>f — k 2 \«m 2 ) will contract to a final vorton state which is strongly 
in the magnetic or electric regimes with \u 2 — k 2 \ ~ 0(m 2 ). Even if the initial 
state was fine-tuned very close to the chiral case cuf ~ kf, then the final state 
always will be less chiral. 

The last point is to check the required condition that fj 2 — > when 

uo 2 — k 2 ~ 0. But in this case, [i and £4 tend to their chiral limits, and we 
have already proved in eq. (26) that in this limit, this quantity is positive, 
£ k=UJ=Q > . So this in turn ensures that the condition for the existence of a 
vorton holds for the whole magnetic and electric regime, except in the limiting 
cases where Eqjv — or Sg^r = 00. 

Limiting cases Q = and N = 0: 

Let us first consider the first of these cases, the so-called spring [5,14], with no 
charge on the loop. It is possible to study analytically the behaviour of these 
objects in this model, if one further assumes that there is no backreaction of 
the condensate on the vortex field. This is a legitimate first approximation, if 
we consider the profiles shown on fig. 1. Then, with our ansatze (22), we are 
able to differentiate (28) with respect to k (note that Sg^r = 0), from which 
we obtain a polynomial in v = —k 2 /k 2 : 
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p( \ - i 1 - V ) 2 dS ^ v 3 / , 5A a v x 2 



+ Y S 4o)^-/i+^S 4o . (29) 



No springs will form if P is always negative. Since we are restricted to the 
interval -1 < v < and as both P(0) < and P(-l) < 0, this will be the 
case if the derivative P' never vanishes between and —1. Now P' is a second 
order polynomial, and some algebra shows that it has a negative minimum for 
v > 0. Thus, it is enough to check that P'(0) = 2/i - ^S 4o > (it can be 
seen that P'(— 1) > by direct computation). To prove this, we note that: 

T^T^^ 2 (30) 

where we have assumed that 4f ^ < 2, and that the string energy density 
is not too different from its critical coupling value. This indicates that for 
reasonable parameters, spring formation is excluded. The behaviour of the 
energy functional is shown in fig. 5. Of course, for the gauged case it is possible 
that additional terms from magnetic pressure could stabilise the configuration 
from collapse [5]. However, since these are logarithmic corrections, such gauged 
springs could only exist on astrophysical scales. 

We have also studied the energetic behaviour of the Q-loop, with TV = 0. In 
this case, we can rewrite the energy functional as: 

£q = 0*- ^ 4 )-^ + 2uQ. (31) 

The typical behaviour of this functional as u varies is shown on fig. 5. As 
can be seen, there does not appear to be any stable configuration, since the 
minimum of the energy is for u = lu c , that is for L = 0. 



5 The stability of the superconducting current 



We have seen in the previous section that the model derived from (1) seems 
to lead to stable loop solutions, but the issue of classical stability has not yet 
been addressed in a full analytic treatment (for the quantum stability, refer 
to [15]). This analysis is inspired in part by a heuristic argument in ref. [16] 
suggesting that classical instabilities could develop in currents with sufficiently 
high winding. Because of the Lorentz invariance of the theory, we have only 
3 cases to consider: pure current (u — 0), pure charge (k — 0), and the chiral 
case uj — k. In this section, we wish to study the perturbations around the 
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Figure 5. This figure shows the behaviour of the energy functional for the spring 
(E is plotted as a function of k 2 ) and the Q-loop (E is then a function of lu 2 ). 
It seems that no stable configuration of finite radius can be expected from these 
configurations. 

pure current solution. To proceed, we decompose our field in the following 
way: a = Sk(r)(e lkz + v(z)), where is the radial part of the unperturbed 
solution and v is an arbitrary perturbation. 



The change in S a is easily computed: 

5 2 E a = Jdz (£ 2 |V z z/| 2 + 7M 2 + ^X 4 [e-* fc V + e lkz v + } 2 ^ , (32) 

where 

1 = JdS {\V r S k \ 2 + ^Si + - \KvI)SI) =S a + - k 2 H . (33) 

Now to simplify the expression (33) for 7, we can make use of our analytic 
result for £ a (15) to find that 7 = — /c 2 S. We can now rewrite the eigenvalue 
equation (following the method used in [17]): 

-E 2 V^ + 1V + y U[v + v + e 2tkz ] = \v . (34) 

To simplify this equation further, consider expanding our perturbation v in the 
form v = e lkz [u\ + iu 2 ]. Substitution leads to the system of coupled equations: 

-EV^i + 2EW,« 2 + (A; 2 S + 7 + A ff S 4 - \) Ul = , (35) 

-£V 2 m 2 - 2EfcV*ui + (A; 2 S + 7 - \)u 2 = . (36) 
To solve this system, we expand in Fourier modes, Ui = die tpz , and we use our 
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simple expression for 7 to obtain the following set of linear equations: 

(Sp 2 + A a £ 4 - A p )ai + 2i(k.p)^a 2 = , (37) 



(Sp 2 - X p )a 2 - 2i(k.p)Ea 1 = . (38) 
The vanishing of the determinant of this system yields our eigenvalues: 



\ p = p 2 S + ± 



^S 4 ) +4S 2 (M 2 



(39) 



Instabilities will occur when one of the possible X p is negative, which is equiv- 
alent to the condition 

SV + 4p 2 S(^S 4 -A; 2 S) <0, (40) 

which leads to: 

— £ 4 - A; 2 S < (41) 
4 

With our ansatz (19) for the moments S 2 , £4, we can use this result to de- 
termine the critical value of ki nst above which a current becomes unstable: 

^sflr (42) 

where a = 4f-£ 4o /£ and k\ is given by (9). Typically, the coefficient a is 
of order unity, and so we can expect A; 2 nst ~ k^/2, in good agreement with 
numerical estimates. We contrast this quantitatively with the heuristic result 
given in ref. [16] which suggested that the onset of instability was indicated 
by the condition |^ < 0; our numerical estimates indicate that the precise 
criterion (42) is significantly lower. 

Note also that through (40) , we are able to predict the associated wavenumber 
p and hence the typical lengthscale of the perturbation for a given unstable 
current with k > k inst . Again, this is in good agreement with our simulations. 
(Strictly speaking, a more accurate analysis should take into account the ef- 
fects of charge conservation; we will discuss this more elaborate calculation 
elsewhere [18], but we note that it yields the same result (42).) 

We point out that this analysis has been carried out by looking only at the 
condensate part of the energy. Therefore, it is valid in a broader class of models 
than the one we are considering here, since the vortex is somehow irrelevant. 
In particular, this instability may appear in any model exhibiting a second- 
order U(l) phase transition, and could be applied in e.g. superconductivity 
studies of resistive transition. 
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t = 54 



Figure 6. Modulus (left column) and real part (right column) of the condensate along 
the string axis at various stages of its evolution, showing the current instability. We 
clearly see that the condensate looses quanta of winding in this process 

Our three-dimensional field theory code allows us to follow the development 
and evolution of these perturbations and their consequences for the structure 
of the superconducting string. From fig. 6, we see that the condensate becomes 
pinched as the instability develops; if the winding number is sufficiently high, 
the instability will force the condensate down such that a = at a localised 
point along the string. This allows the string to lose quanta of winding and 
lowers k in this region. The condensate then bounces back and slowly relaxes 
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into the stable current configuration, radiating the energy difference between 
the initial and the final configuration. These numerical results are provided by 
way of illustration but a more extensive study will be reported shortly with 
further details about the simulation code [18]. 



6 Conclusions and prospects 



In this paper we have studied the behaviour and effects of the currents and 
charges on superconducting strings and loops. We have derived exact formulae 
for the energy of a loop, which exhibits a generic divergence at u 2 — k 2 = 0, 
therefore proving that the chiral case is not an attractor, but rather a repeller. 

Since for stable configurations we typically have E ~ ^qn = Q/N, from 
(27) we can expect the final size of the loop to be very much smaller than 
its original size (for realistic initial conditions). Of course, we have not taken 
into account the vortex-antivortex interactions on scales small compared to 
the string width, so the inevitability of vorton formation is subject to this 
important caveat. In addition we conclude that springs and Q-loops are not 
allowed in this theory. 

Studying the stability of the current on a straight string, we have also seen 
that the superconducting regime is unstable when the winding is too high. 
By Lorentz-invariance, it is easy to see that this will be the case if m 2 v = 
k 2 — oj 2 > kf nst7 where ml ~ k 2 ~ u 2 . (We believe that the gauged case will 
exhibit the same generic features as these global currents, though with the 
small quantitative differences already discussed.) 

Now, from our analysis of the vorton state, we know that equilibrium typically 
will be achieved when £ ~ Eqat. Using our ansatz (22), this is equivalent to 

v - v 77 • ( 43 ) 

To ensure stability, we have to impose v > —k 2 nst /m 2 which, with (43) and 
some algebra, leads to the condition 

fry) 2 L>2 

^qn > \ - rjF S c . (44) 

m o K inst 

Loops that do not satisfy (44) will be unstable and lose quanta of winding. 
Hence, Sq^t will increase, and the loop may reach a stable state. This process 
is associated with energy radiation, which would be interesting to quantify to 
determine possible observational signatures of this phenomenon. 
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Finally, the perturbative stability analysis carried out for superconducting 
currents in the magnetic regime has been extended to the chiral and electric 
cases, and appears to establish their stability. However, this issue is more 
subtle, and the analysis will be published separately [18]. 
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